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ABSTRACT 

Like the majority of spiral galaxies, NGC 6155 exhibits an exponential surface brightness profile 
that steepens significantly at large radii. Using the VIRUS-P IFU spectrograph, we have gathered 
spatially resolved spectra of the system. Modifying the GANDALF spectral fitting routine for use 
on the complex stellar populations found in spirals, we find that the average stellar ages increase 
significantly beyond the profile break radius. This result is in good agreement with recent simulations 
that predict the outskirts of disk galaxies are populated through stellar migration. With the ability 
to bin multiple fibers, we are able to measure stellar population ages down to ^lv '^ 24 mag/sq arcsec. 
Subject headings: galaxies: formation — galaxies: structure — galaxies: stellar content — galaxies: 
individual(NGC 6155) — galaxies: spiral 
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1. INTRODUCTION 

The nature of how and where stellar disks end can shed 
light on the limits of star formation and the importance 
of secular evolution. Stell ar disks tend to be embedded 
in niuch larger HI disks (jBosmal 119811 : iBroeils fc Rhed 
119971 : [Begum et a"Lll2005[ ). showing that while there are 
baryons in the outskirts of galaxies, they are not being 
converted into stars efficiently. It is often assumed that 
there is a surfac e density threshold for star formation 
(|Kennicuttill989D which would provide a natural mecha- 
nism for truncating the stellar component of galaxies. 

The surface brightness profile s of spiral gala xies are 
commonly fit with exponentials (|Freemanlll97Clt ). While 
the exponential function is a convenient description, 
few astronomers would suggest that stellar disks ex- 
tend to infinite radii. Observations of edge-on sys- 
tems originally suggested that stellar disks have well 
defined truncation radi i (Ivan der Kruit fc Searld 119811 : 
iKregel fc van der Kruitl 120041) Using large s amples 
of nearly face-on galaxies, iPohlen fc Trujilld (|2006f l 
found a variety of beh aviors at large galactic radii. 
IPohlen fc Truiillol ()2006[ ) find that only 10% of disks are 
well fit by a single exponential light profile, while in 60% 
of their sample the surface brightness profile is well fit 
by a downbending broken exponental, with a break ra- 
dius 1.5-4.5 times the inner scale length. The final 30% 
of their sample are described with a broken upbending 
exponential profile. Unlike the edge-on studies, they do 
n ot find evidence for sharp truncations. 

iRoskar et al.l (|2008D present a possible formation mech- 
anism for surface brightness profile breaks. They find 
that a star formation threshold seeds the onset of a 
truncation, but radial migration populates the disk be- 
yond the break producing a down-bending broken pro- 
file. Since the stellar migration is a random walk, 
only ancient stars have had enough time to get to 
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very large radii beyond where star formation shuts 
off. These simulations thus predict that the break ra- 
dius should correspond with an increase in the average 
stellar ages. The formation and appearance of outer 
disks has since been investigated in cosmological sim- 
ulations by several group s who also find an upturn in 
age at the break ra dius fS anchez-Blazguez et al.l 120091 : 
iMartinez-Serrano et al. 2009il. However, the relative im- 
portance of radial migration and in-situ star formation 
in a fully-cosmological setting remains an open question. 

While simulations find that profile breaks can be 
caused by a combination of a star formation threshold 
and stellar migration, observations suggest that star for- 
mation can pro ceed even in th e outskirts of disks. Deep 
Ha imaging by iFerguson et al.l (Il998b.) s how s igns of ex- 
tended star formation. iThilker et al.l (|2007l ) also find 
that ~ 30% of spiral galaxies have UV emission extend- 
ing beyond thei r traditional star fo r matio n threshold ra- 
dius. Similarlv. iHerbert-Fort et al.l (|2009|) find clustering 
of point sources around NGC 3184, implying star cluster 
formation in the outer disk. 

Recent observations of resolved stars have found 
a variety of r esults for stellar truncation regions. 
I Williams et~all (|2009al ) find that in low-mass M33 the 
the age gradient changes from negative to positive a t the 
truncation radius. Meanwhile, Ide Jong et al.l ([2007' ) find 
that the truncati on region in NGC 42 44 is the same for 
stars of all ages. Ide Jong et al.l (|2007D interpret this as 
a sign the truncation is the result of a dynamical inter- 
action, but these results were subsequently shown to be 
consistent with the radial migration model ([Roskar et al.l 
1 20081 . See iVlaiid pOlO ) for a review of the recent ad- 
vances in studying the outskirts of spiral disks. The 
observational results point to an interesting conundrum 
that while many studies find active star formation in the 
outskirts of disk galaxies, the stars in the region are often 
older than the inner disk. 

To date, truncation regions have been studied with 
broadband colors and resolved stars. There are a lim- 
ited number of systems we can resolve, and faint broad- 
band surface photometry is dominated by systematic er- 
rors associated with flat-fielding and sky subtraction. 
While there have been searches targeting emission lines 
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(jFerguson et al.lllQQSaHChristlein et al.ll2010[ ). in this let- 
ter we present the first spectroscopic study of stellar con- 
tinuum spectral features across a surface brightness pro- 
file break. 

2. OBSERVATIONS 

Observations of NGC 6155 (RA 16:26:08, DEC 
-h48:22:01, J2000) were taken with the Visible Integral 
fiel d Replicable Un it Spectrograph Prototype (VIRUS- 
P) (|Hill et alJl2"008l ) on the McDonald Observatory 2.7m 
Harlan J. Smith telescope on April 9 and 11, 2008. 
VIRUS-P is a 1.7'x 1.7' field-of-view IFU spectrograph 
with 246 4.3" diameter fibers. Observations were taken 
in a three-step dither pattern to attain full spatial cover- 
age. We observed four 20 minute exposures at each dither 
position for a total of four hours of exposure time. The 
data were reduced with the VACCINE package (Adams 
et al. in prep) in combination with custom IDL rou- 
tines. Wavelength calibration was performed using ob- 
servations of Hg and Cd arc lamps. Twilight sky ex- 
posures were used to generate flat-field frames. We ob- 
served the stars Feige 34 and BD-l-33 2642 to establish 
flux calibration. The final result is a data cube with 738 
spectra filling a 1.7'x 1.7' field of view and a FWHM 
resolution of 5.3 A with 2.2 A pixels and a wavelength 
range of 3550-5580 A. We mask a 10 A region around the 
bright 5577 A O I skyline. 

Using the T^-band surface brightness measured from 
our data cube we find an inclination of 46 degrees and 
position angle of 135 degrees (East of North). These val- 
ues are consistent with the measured velocity field. The 
radial surface brightness profile is well fit with a broken 
exponential with a break radius of 36", where the scale 
length changes from 11.6" to 7.7". We fit the region of 
5"< R <60" to avoid any centra l bulge and outer regions 
of low SNR. IPohlen fc Tr^millol 12006) fit the SDSS im- 
age of NGC 6155 and find a break radius of 34" and scale 
lengths of 12.2 (13.1) and 8.1 (7.3) in t he r (p) band. Our 
values are an excellent match to the IPohlen fc Truiillol 
([2006) fits, especially considering our lower spatial reso- 
lution. Figure [1] shows our fiber photometry, the velocity 
field measured from gas emission lines (iif3, [Oiii], and 
[On]), and our radial profile fit. 

We re-binned the spectra to a common velocity. Spec- 
tra were then co-added in 4" elliptical annuli of constant 
galactic radius to reach adequate SNR at large radii. 
Fibers that contained bright stars or background galaxies 
were masked. Our final binned spectra have SNR/A~320 
in the inner regions and ~ 45 beyond the break radius. 
We assume a Virgo infall corrected redshift to NGC 6155 
of 2549kms~i(38.3Mpc). 

3. STELLAR POPULATION SYNTHESIS 

Historically, prominent individual spectral absorp- 
tion features have been used to constrain the flux- 
weighted average age of stellar populations. Systems 
based on stellar absorption features like Lick indices 
have proven useful for old elliptical stellar po pulations 
(e.g.. IWorthev et al.|[l99l iTrageret "all [2000 ). but be- 
come difficult to interpret when there are multiple stel- 
lar populations or e xtended star formation histories 
(|Serra fc Traged[2007l) . 

There are now a host of packages developed to 
measure star formation histories (SFHs) from integrated 



spect ra (e.g.. lOcvirk et al.l 120061: iMacArthur et al.l 
200l iToieiro et all 1 20071: iKoleva et all [200a 
Chilingarian et al.l I2007D . All of these codes at- 
tempt to model observed galactic spectra as a linear 
combination of single-burst simple stellar populations 
(SSPs). While it should be possible in theory to recon- 
struct a star formation history as a series of SSPs, in 
practice, we find that the problem becomes degenerate 
at low SNR. There is also a strong degeneracy between 
including dust extinction and the presence of an old red 
stellar population. In general, we find that it is very 
difficult to constrain an average stellar age using a linear 
combination of SSPs. 

Unlike spheroidal systems that can be well described as 
having only old stars, spiral disks overwhelmingly display 
a mix of very old and newly formed stars. Observations of 
resolved stars in nearby spiral galaxies show that st ar for- 
mation histories are extended at all radii (Williams et al.l 
I2009allbl : iGogarten et al.ll2010HBarker et al.ll2007D . Sim- 
ilarly, the solar neighborhood in the Milky Way contains 
young stars, intermediate age stars like the sun, and very 
old thick disk stars. 

Motivated by the observational results that even the 
outskirts of galaxies have a mix of stellar ages, we 
have constructed a large library of extended SFHs and 
find the best fitting single spectrum. We use the lat- 
est stellar synt hesis package from Chariot & Bruzual 
(jCharlot fc Bru zual 2007; Bruzual 2007) to generate a 
suite of composite stellar population spectra. For all the 
models, we assume an exponential star formation history 
of the form 



V'(t, [Z/H]) = e-*/^M0yr- 



(1) 



where "0 is the star formation rate and r is the star form- 
ing time-scale. We include both positive and negative 
values of r. All the models are computed at an age of 13 
Gyr. We generate 1664 total model spectra, with values 
of r ranging from -200 Gyr to -0.25 Gyr, and 0.5 Gyr 
to 200 Gyr. Negative values of t represent an increasing 
star formation rate. Small absolute values of r repre- 
sent extreme stellar populations (dominated by very old 
stars or newly formed stars), while large absolute values 
of T are effectively constant star formation histories. Our 
models span mass-weighted ages of 15 Myr to 12.7 Gyr. 
The metallicity values run from Z=10~'* to 0.1 (with so- 
lar defined as Z=0.02). We use the Padova 1994 steUar 
evolutionary tracks and a Chabrier initial mass function 
(IMF). 

We do not include any metallicity evolution in the 
models, thus our fits only return a flux-weighted aver- 
age metallicity. The advantage of our method is that it 
runs fairly quickly because of the limited parameter space 
and the SFHs are motivated by observations of spatially 
resolved stellar populations. The disadvantage is that 
our metallicity measure is fairly crude since we do not 
include chemical evolution. Also, we can not expect our 
flts to be robust if a galaxy has undergone strong episodic 
bursts of star formation. We do not include a-element 
variations in our spectral models, however as a late-type 
spiral, we would not expect NGC 6155 to be particularly 
a— enhanced. We emphasize that our fitted SFHs are a 
measure of the current stellar populations, after radial 
mixing and azimuthal binning, and not a measure of the 
actual star formation history at a fixed radius. 
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Fig. 1. — Upper Left: SDSS image of NGC 6155 (North is up and East is right). Upper Right: Surface brightness measured from 
VIRUS-P fiber spectra. Points are the same size as the VIRUS-P fibers. Fibers containing stars or bright background galaxies have been 
masked. Lower Left: Velocity field measured from emission lines. Lower Right: Surface brightness profile as measured by individual fibers. 
A green line shows the best-fit broken exponential. The dashed green line shows the result of extrapolating the inner and outer region fits. 
A yellow curve shows a running median with an 5" window. The dashed line shows the level where the SNR of a single fiber is unity. In 
each panel, a red curve marks the best-fit break radius. 
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We use the Ga s AND Absorptio n Line Fitting (GAN- 
DALF) package (jSarzi et al.ll2006D along with the penal- 
ized pixel fitting (pPXF) ("Capp ellari fc Emsellemlhool) 
code to find the best-fitting model for each binned spec- 
trum. The pPXF code fits a stellar velocity and disper- 
sion which is then held constant for the rest of the fit. 
GANDALF multiplies model spectra by a low-order poly- 
nomial (to correct for dust extinction and flux calibration 
mis- matches), and simultaneously fits specified emission 
lines as Gaussians. Figure [2] shows an example of our fit- 
ting procedure. The use of the polynomial ensures that 
the fit matches the spectral absorption features, and the 
overall tilt or broadband color does not affect the fit. 
Because dust extinction has a negl igible effect on th e rel- 
ative depth of absorption features (|MacArthuHl2005 !) , our 
fitting method does not require correcting the spectra for 
internal or foreground extinction. 

The uncertainty in our fits should be dominated by the 
systematic difference between the adopted smooth SFH 
and the bursty SFHs of real galaxies. To test the errors 
associated with SFH mismatch, we have generated inte- 
grated spectra using the SFHs measured from resolved 
populations in NGC 300 (Gogartcn ct al. 2010) and fit 
them with our r— SFHs. The inner regions of NGC 300 
have SFHs strongly peaked at an age of 12 Gyr while 
the outer regions have constant or increasing SFRs. The 
radial bins of NGC 300 have average ages ranging from 
5.2-9.6 Gyr. Fitting our exponential SFHs, we are able 
to match the relative ages of each radial bin in NGC 300 
and have an absolute age dispersion of a = 0.6 Gyr. The 
largest source of error seems to be recent bursts of star 
formation causing our fitting routine to converge on ages 
that are younger than the actual mass-weighted stellar 
population. 

Figure [2] shows an example of our data and fitting pro- 
cedure. The fits become noticeably worse blueward of 
^4300 A. There are several possible reasons for the poor 
fit in the blue. The instrument throughput in the blue 
during these observations was low, resulting in poorer 
SNR. The templates may need improvement in the blue. 
Finally, the poor fit could also be a systematic error left 
over from imperfect sky subtraction. We have experi- 
mented with masking this region and find the fits remain 
basically unchanged. We emphasize that our errors are 
clearly non-Gaussian and are dominated by the system- 
atics of template mismatch and SFH mismatch. There- 
fore, rather than adopt the standard error bars derived 
from 8^ — 2.3 contour, we plot the larger error bars 
derived from our experiment fitting the observed SFH 
of NGC 300. We also note that a recent burst of star 
formation could bias our measurements to return values 
that are too young. 

4. RESULTS 

The results of our stellar population fits are plotted in 
Figure [3] The central region is best fit by a decreasing 
star formation rate, with a mass- weighted age of ^ 7 Gyr. 
The region at a radius of 15" (^2.8 kpc) is dominated 
by very young stars and active star formation. This is 
clearly visible in the SDSS image as a region dominated 
by a star forming ring or inner spiral arms. Beyond the 
star forming region, the stars remain young and domi- 
nated by recent formation. Past the surface brightness 
break radius at 34" (~ 7 kpc), there is a sudden change. 



and the stars are best described with older stellar popu- 
lations and declining SFHs. This is an impressive confir- 
mation of the hypothesis that stars beyond the break ra- 
dius are populated by stellar migration and should show 
an upturn in stellar age. 

Our fits return a metallicity [Z/H]~ —0.2 for most of 
the region before the profile break. Beyond the profile 
break, the metallicity drops to [Z/H]~ —0.5. While this 
is consistent with the stellar migration hypothesis, we 
caution that this is still only a crude measure of the stel- 
lar metallicities. 

5. DISCUSSION 



IPohlen fc Tr^lililol ([20061) describe NGC 6155 as lop- 
sided and classify the galaxy as a Type II-AB (a down- 
bending profile in an asymmetric disk). While our fiber- 
photometry shows slight lopsidedness, the velocity field 
is quite regular and there are no nearby neighbors. All 
the galaxies within a 1 Mpc projected radius and similar 
redshift of NGC 6155 are at least 3 magnitudes fainter. 
It therefore seems unlikely that the profile break in NGC 
6155 is the result of a dynamical interaction with a neigh- 
boring galaxy. 

Our result that average stellar ages increas e beyond the 
break radius is consistent with the results of lBakos et al.l 
(|2008D who stacked multiple SDSS images and found the 
g — r color became redder beyond the profile break. This 
is also similar to the age gradient observed in M33, with 
the age decreasing out to the profile break where the gra- 
dient reverses and the stellar population b ecomes older 
beyon d the break ([Williams et al.l l2009a|) . iBakos et al.l 
([2008() claim that while their stacked surface brightness 
profiles shows a break, the stellar mass profile is de- 
scribed by a single exponential. Our model mass-to-light 
values are not well constrained beyond the profile break, 
however, the increase in age we detect does suggest the 
M/L is larger in this region. Therefore we expect the 
stellar mass profile to be smoother than the light profile. 

The increase in stellar age we measure is very similar 
to the simulations of Roskar ct al. (2008) where the av- 
erage stellar age jumps from 4 to 6 Gyr. The best fits 
beyond the break radius have r ^ — 9.5 Gyr, consistent 
with significant ongoing star formation in the region. We 
therefore take our results as evidence that while stellar 
migration can populate outer disks, our fits are consis- 
tent with some low-level in situ star formation as well. 

Along with the change in stellar age, the surface bright- 
ness break also corresponds with a lack of line emission 
(Figure [TJ lower left). This could be a simple coinci- 
dence where the break radius is also where the emission 
becomes faint enough that we no longer detect it, but it 
is also consistent with the theory that the break radius 
represents the location of a star formation threshold. 

It should be pointed out that our observations are also 
consistent with a radial change in the stellar IMF. If 
the outer disk only hosts lower-mass molecular clouds, 
we could expect star formati on to be unable to pro- 
duce the h ighest mass stars ([Kroupa fc Weidneil 120031 : 
iKoppen eTa l. 2007). Because young massive stars con- 
tribute heavily to the overall luminosity, their absence 
would result in a young stellar population that spec- 
troscopically appears older. A change in the IMF can 
also explain the lack of emission lines beyond the break 
([Pfiamm-Ahcnburg fc Kroupa, 2008 ). 
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Fig. 2. — An example of our stellar population fitting procedure. Top: The computed x surface. Each point represents a unique spectral 
model. The minimum along with contour levels of Ax =20 and 200 are plotted. Bottom: The observed spectrum (black) with the laest-fit 
model (red). Residuals are plotted. 
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In a forthcoming paper, we will present a larger sam- 
ple of similar observations. Our larger sample includes 
galaxies with no surface brightness profile break as well 
as up-bending profiles. We also plan to further refine our 
spectral synthesis template fitting method to improve the 
metallicity measurements of the stellar populations. 
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Fig. 3. — Top: Radial mass-weighted age profile for NGC 6155. Red circles are used for decreasing star formation histories and blue 
triangles are used for increasing histories. The surface brightness break radius is marked with a red vertical line. Error bars represent our 
estimated systematic uncertainties. Middle: The traditional spectroscopic age indicator D4000. Bottom: Emission line corrected Balmer 
absorption equivalent widths. We have not matched resolution to formally move onto the Lick system and show the Blamer strengths 
for illustration. The D4000 index as well as the Balmer lines all show changes at the profile break consistent with the stellar population 
becoming older. The bright star formation ring around R=15" is also clearly measured as younger. 



